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The framework of the expansion by subgraphs is used to compute asymp¬ 
totic expansions for the vector and the tensor currents in the limit of large 
quark masses. We use the results to obtain an estimate for the influence of 
heavy quarks on the nucleon electromagnetic and tensor form factors. 


1. Introduction 


The intrinsic heavy quark content of light hadrons is a fundamental property of QCD 
[1, 2, 3]. The influence of intrinsic heavy quark content on the structure of the nucleon 
can be investigated by several methods. Much information about the nucleon structure 
like the charge distribution and the anomalous magnetic moment can be obtained from 
form factors. These form factors parametrize the expectation values of bilinear fermion 
operators in a single-nucleon state. To describe the influence of intrinsic heavy quarks 
on the form factors one should therefore consider those bilinear fermion operators that 
contain heavy quark fields, it is natural to employ a large mass expansion for the in¬ 
vestigation of such operators. One way to perform large mass expansions is the heavy 
quark mass expansion (HQME) which was used in [4] and [5] to obtain an expansion of 
the vector and the tensor current of heavy quarks in terms of gluon field operators. 
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One observation that the HQME-approach is not complete was made in [6]. This ob¬ 
servation concerns the influence of muons on the anomalous magnetic moment /xe of the 
electron in QED. If the muon mass is treated as large parameter (compared to the elec¬ 
tron mass) this influence can be described by using the HQME to expand the muonic 
vector current in terms of photon held operators. The first photon field operators in the 
expansion that contribute to fxe are of order l/mf^ [7]. However, from the evaluation of 
the corresponding Feynman diagrams one can show [8] that the contributions of muons 
to He start at order 1/m^. Due to the similarity of the calculations, the same problems 
can be expected for the HQME of the vector current of heavy quarks in QCD. In [6] 
it was argued that the problem in QED arises because the 1/m^-suppressed photon op¬ 
erator has to be renormalized by 1/m^-suppressed electron operators. That settles the 
problem in this particular case, but it would be desirable to have a systematic procedure 
to find those additional operators that complete the operator product expansions of the 
HQME. 

To give such a systematic procedure, we use the language of Feynman diagrams instead 
of that of path integrals that was used to derive the HQME results. In this language, 
rigorous results about the influence of virtual heavy particles are known for a long time. 
In Ref. [9] the decoupling theorem was proved: one can find renormalization schemes, 
in which, to leading order in the mass of the heavy particles, the effects of the heavy 
degrees of freedom can be neglected altogether. Afterwards much work was devoted to 
the problem of how to obtain a systematic expansion that allows to take into account 
the effects of heavy particles both on the diagrammatic level and on the level of the 
Green functions (see for example Ref. [10], Ref. [H])- The first treatments were 
done within a momentum subtraction scheme because here the close connection between 
asymptotic expansions and renormalization was most apparent. However, due to the 
fact that dimensional regularization is nowadays the preferred method for regularization 
both because of its technical simplicity and its preservation of gauge invariance, it was 
desirable to derive similar results within the minimal subtraction scheme. In the MS- 
scheme decoupling is not straightforward: not all effects of heavy particles are suppressed 
by powers of their masses. Instead, one has to absorb effects of heavy particles in the 
physical coupling constants of the theory [12]. Subsequently a systematic procedure 
for the asymptotic expansions within the MS-scheme was developed (see for example 
[13]). Finally the technique of ‘expansion by subgraphs’ [14] was developed, for which it 
was proved that one obtains an expansion that is free of infrared divergencies and that 
has the correct 1/M behaviour. This technique works on the diagramatic level in the 
first place. However, with the help of the counterterm technique these results could be 
immediately generalized to an operator product expansion. 

Here we will use this technique to analyze the validity of the HQME-approach. This is 
done for the especially important examples of the vector current and the tensor current. 
The Section explains the basic idea of asymptotic expansions. In the following section 
the expansion is carried out for the vector current and the tensor current. Finally, we 
use the expansion of the vector current to estimate the influence of charm quarks on the 
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e.m. and tensor form factors of the nucleon. The reader interested in the final result can 
directly go to Section 5. Applications for the intrinsic charm in nucleon form factors are 
discussed in Sections 6 and 7. 


2. Asymptotic expansions in heavy quark masses 


Perhaps the simplest way to illustrate the general idea is to cite the power counting 
theorem for the dependence of a Feynman diagram on the mass of a heavy particle from 
Ref. [15]: 

Let r be a diagram with I loops. The asymptotically irreducible diagrams 
(Al-diagrams) of r are those connected subdiagrams that contain heavy lines 
and cannot be made disconnected by cutting a single light line. Then the 
asymptotic behaviour of r is bounded by 

log(m)* 

where S is that set of disconnected AI graphs (spinney) containing all heavy 
lines which has the highest degree of divergence. 

From this it is immediately clear that by introducing an operator that lowers the degree 
of divergence of the Al-diagrams sufficiently (like the R-operation in momentum sub¬ 
traction schemes lowers the degree of divergence to make IPI graphs convergent) one 
can generate a remainder that is suppressed by an arbitrarily chosen power of the heavy 
mass. The difference between this remainder and the original (renormalized) diagram 
will then constitute the correct asymptotic expansion for the diagram. 

This procedure is systematically described in [13]. We only give the main result. 

We first have introduce some notation. We denote the large mass in which the expansion 
is done by M, the external momenta are collectively denoted by {p}. Let r be an arbitrary 
graph and let the diagram 7r({p}, m,M) be the corresponding analytic expression. The 
degree of divergence of the diagram is w-y. Let be the set of all Al-spinneys 

(set of mutually disjoint AI subgraphs) and I^/s the diagram /p with all elements of 
the spinney S shrunk to a point. Ip/s ° 11765^7 is the diagram in which the vertex Fy 
replaces the subgraph 7 for each 7 g S. R"”/r/5 is the R-operation that only acts on 
those parts of the diagram that do not contain vertices to which the graphs of S were 
shrunk. Finally, My'*' performs a Taylor expansion of the diagram / 7 ({p},m,M) in the 
momenta {p} external to the graph 7 and the light masses m up to the order ay. 

In [13] it was shown that for a-y = a + oj-y one has, up to O (j-terms, 

RIr{{p},m,M) it""7r/5({p},m) o [Q My^ i?/y({p}, {<j},m,M). (1) 

seSyij(r),s/o j&s 

It was also proven that this result is free of artificial UV- and IR-divergences, that is, 
the divergences of this expression are the same as those for the original diagram. To 
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summarize; To perform an asymptotic expansion in large masses on has to find all AI- 
spinneys of a diagram r. In each subgraph 7 of this spinney, one performs a Taylor 
expansion in external momenta and light masses. The result is a new vertex factor 
Vy that is reinserted in the original diagram to replace the subdiagram 7 . Then the 
unaffected parts of r are renormalized as usual. 

The transition from the diagrammatic to the operator level is now fairly obvious: If one 
computes the matrix element of an operator containing heavy particles one will obtain 
Feynman diagrams containing Al-subgraphs. The expansion of these diagrams in terms 
of the heavy mass amounts to replacing the AI subgraphs by their Taylor expansions in 
light masses and external momenta. The terms of this expansion serve as new effective 
vertices in the complete diagram. The complete dependence on the heavy mass is in 
these effective vertices. Thus, if we reexpress the effective vertices as vertices due to 
local operators in the light degrees of freedom, all matrix elements of the heavy operator 
can be described by matrix elements of light operators. 

Now, if one wants to expand the heavy operator O of dimension dim{0) to the order o, 
on the diagrammatic level one has to expand each Al-subgraph to the order = cuy + a. 
Since in theories with dimensionless coupling constant ujy = dim{0) — Eb — ^Ep {Ep, 
Ep denote the number of external boson- resp. fermion-propagators), it is enough to 
consider diagrams with 0 < a - 1 - dim{0) - Ep - ^Ep. For example, at the order 1/M^ for 
dimip) = 3 only diagrams with Ep < 5 and Ep < 2 (remember that external fermion lines 
have to come in pairs) have to be taken into account. In gauge theories. Ward-identities 
will lower the actual degree of divergence for diagrams with external vector bosons. This 
reduces the number of diagrams that have to be considered. 

In this paper the vector current and the tensor current are considered. For these opera¬ 
tors an OPE is done to the order 1/M^ for the vector current and to the order 1/M for 
the tensor current and up to the order in the strong coupling constant. The reason 
for the order is the fact that this is the lowest order with purely fermionic operators. 


3. The Al-diagrams for fermion-bilinears 


We first give a general constraint for the structure of the expansion (see the appendix 
of Ref. [10]). For the argument that follows we use a regulator that admits the usual 
definition of 75 . Since the heavy particles should not occur as external particles, every 
Al-diagram will contain heavy fermions only inside closed loops. If we switch the sign 
of M, we have 


Sp{p,-M)= -j^Sp{p,M)'r^. 
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Now, the factors of 7 ^ can be absorbed in the vertex factors via 


5 u 5 
_ "Y ■“ 

5 uu 5 
-7 0-^7 = 


This cancels all 7 ® in the loop, leaving the diagrams with vector current insertion as 
they are and changing the sign of the diagrams with tensor current insertion. One can 
conclude that in the expansion of the vector current only even powers of M can occur 
while the expansion of the tensor current contains only odd powers of M. 

We can further restrict the types of relevant diagrams. Only gluons can couple directly 
to the heavy quark loop. This means that each AI diagram will have a subdiagram that 
can be considered as off-shell matrix element of 


{x,pi ,... ,p„) := {0| TO{x)Af^^ (pi)... (p„) |0). 

Matrix elements with one external gluon are zero by color conservation. Let us consider 
the case of two external gluons. Under charge conjugation, we have 



Thus, using global color conservation results in 


{x,pi,P2) = {0|T(!i(a;)A^i(pi)T7(p2)|0) 

= ^ (0| TO{x)trc (pi) (p 2 )) |0) 

{X,P 1 ,P 2 ) (-1)^ {0|rc>(a;)trc (^(^^i)^(pi) (A'' 27 (p 2 )^ |0) 

= - ^ (0| TO{x)trc (Pi) (p 2 )) |0) = 0. 

This shows that for the vector current and the tensor current one needs to consider only 
diagrams where the fermion loop is connected to the rest of the diagram by at least three 
gluon propagators. 

One can also find the only possible color structure for the matrix elements with three 
external gluons: From global color conservation we see that there are only two relevant 
color structures for the matrix element of the current with three external gluons (8 (g> 
8 (g> 8 contains two different singlets). These can be taken as the two independent trace 
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structures. In combination with charge conjugation invariance this leads to: 


^ abc 


ix,Pl,P2,P3) = 


(0| TO{x)A'^^ (pi) <2 (pg) |0^ 

- ^ { 0 \TO{x)trc (pi) [yl^2 (P2) (pg)]) |0) 

+ ^ ( 0 | TO{x)trc (pi) {^^^2 (p2), (p3) }) |o) 

^ jabc ^0|rc>(x)trc (pi) [^^2 [p 2 ),A^^ (P 3 )]) |0) 

+ ^ (0| ro(x)trc (Pi) {^^^2 (P2), (P3) }) |0) 

{x,pi,P2,P3) ■ 


Thus, 


(0|TO(x)tre(^'^l(Pl)KMP2),^'^3(P3)])|0> = 0. 

For the case of the vector current one can conclude from QED current conservation that 
there is an additional factor of the momentum entering the diagram at the operator 
insertion. By dimensional analysis, the dimension in the heavy mass is lowered by by 1. 
We will show in section 4.1.1 that if three gluons couple directly to the fermion loop each 
of the external gluons comes with a factor of external momentum. Taking into account 
both kinds of momentum factors, we have the power in the heavy mass ZJ^ ^ ojj - Eg — S 
where Eg is the number of external gluons coupled directly to the heavy loop, <5 = 1 for 
the vector current and <5 = 0 for the tensor current if no more than three gluons couple 
directly to the heavy loop. If there are more than three gluons coupled to the loop, we 
have — S. 

Now, we list the Al-diagrams relevant for us as shown in figure 1. We start with one 
loop diagrams. Following the arguments given above, the diagrams with three external 
gluons [la] give the first contributions for the vector current and the tensor current. 
For the vector current these diagrams have = —4, while for the tensor current they 
have cU^ = -3. Actually these diagrams can be excluded from the computations at the 
order we consider. They are nevertheless, because we want to compare the results to the 
HQME-results given in [5]. The other group of diagrams that is relevant at the order 
considered in our expansion is that with four external legs [lb]. By power counting, 
these diagrams have = —2 for the vector current and = —1 for the tensor current. 
However, our explicit computation and the results of section 4.1.1 showed that there are 
no 1 /m2 or 1/M contributions. This is consistent with the HQME results which predict 
that all one loop contributions are suppressed up to the order 1 /M"* for the vector current 
and 1/M^ for the tensor current. 

Now, consider two-loop diagrams. Since at least three gluon lines must couple to the 
fermion loop, the Al-diagrams that are relevant for us are given by the figures Ic, Id, le 
and If at order pf and by Ig Ih, li and Ij at order p®. 

The three loop diagrams are the first which contain external external fermions. In the 
computation we will only consider one light quark flavor since all light flavors will, up to 
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(a) 1 loop: type I (b) 1 loop: type II (c) 2 loop: type I (d) 2 loop: type II 



(i) 2 loop: type VII (j) 2 loop: type (k) 3 loop: type I 
VIII 

Figure 1: Al-diagrams for the HQME of vector current and tensor current: The blob 
denotes the operator insertion, the double lines denote heavy quarks 


trivial modifications, yield the same result. If only two of the gluons are coupled to the 
fermion line, we have at least = —3 for the vector current and = —2 for the tensor 
current. Thus, the only possibility is to couple all three gluon legs to the fermion line, 
see Fig. Ik. For external ghost fields charge conjugation invariance shows the absence 
of contributions with two external ghosts (due to the ghost number conservation the 
number of ghosts fields must be even): 

G 2 ^ {x,Pi,P 2 ) = (0| TC>(a;)ca (pi) Cb (p 2 ) |0) 

= ^ (0|TC)(a;)cd(pi)cd(p2) |0) 

^ (-l)^(0|TO(x)crf(pi)cd(p2)|0> 

= 0 . 

All nonzero diagrams with more external ghost fields are of higher order than 1/M^ in 
the expansion. 

Note that to all diagrams given in Fig. 1 one also has to compute the crossed diagrams. 


There is one additional restriction for the purely gluonic operators in the expansion: 
The expansion should respect local gauge invariance. Thus, it should contain only field 
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strength tensors and covariant derivatives of them. Since the behaviour of the field 
strength tensor under C is the same as that for the vector potential, the arguments 
given above show that the result for both operators must have the schematic structure 
trc{{Fi,F 2 }F 3 ). For the vector current, current conservation requires that the result 
must be the derivative of some antisymmetric tensor, giving d trc{{Fi,F2}F3). From 
power counting it is now clear that the gluonic contributions to the expansion should 
be of order 1/M^ for the vector current and 1/M^ for the tensor current to all orders of 
perturbation theory. Our computations up to 0{af) confirm these results. 


4. Computation of the Al-diagrams 


We can now use the relation (1) to compute the large mass expansion for the diagrams 
shown in the previous section. There is one additional simplification. For the expansion 
considered here single diagrams may have divergences, but if one takes the sum over all 
diagrams that contribute at a given order in the heavy mass and the coupling constant, 
the divergences cancel. Actually even the sums of all diagrams of a given type (as 
indicated in Fig. 1) are all finite. Thus, it is not necessary to perform the R-operation 
for them. 

To obtain the desired expansion, each Al-diagram has to be expanded in its small masses 
and external momenta. However, it is more convenient to rescale the integration mo¬ 
menta according to I ^ Ml (note that this makes the integration variables dimensionless). 
Since the integrand I{M,{l},{m},{p}) is a homogeneous function (here {R labels the set 
of integration momenta and {p} external momenta while {m} is the set of light masses), 
we have for an n-loop integral 

I{M, {Ml}, {m}, {p}) = { 1 }, [m/M], {p/M}). 

This makes it obvious that we can also expand in 1/M. 

After the Taylor expansion of the of the AI diagrams we obtain several thousand of dif¬ 
ferent tensor integrals and after tensor reduction there are thousands of scalar integrals. 
In the following we explain how the tensor reduction was done and how each type of 
diagram was computed. The computation was done by a Mathematica program which 
made use of the package ‘FeynCalc’ [16]. 

In most cases the expansion of the Al-diagrams given above does not result in scalar 
integrals. This is either because of uncontracted Lorentz-indices in the diagram or be¬ 
cause of scalar products of external and loop momenta occurring in the numerator in the 
1/M-expansion of the integrand. We will explain how the tensor integrals were reduced 
to scalar ones. Since the expansion is in all external momenta, the (tensor) integrals 
that we obtain effectively correspond to diagrams without external momenta; they are 
vacuum diagrams. For these, the tensor decomposition is easy: Only the metric tensor 



is available to represent the tensor structure of a diagram. Thus, an arbitrary integral 
- (all integrals with an odd number of indices vanish) occurring in the expansion 
can be written as 


7 ?Mi -A‘ 2 n = Aig^^ ^^2 . gf^ 2 n-i ^^ 2 n _|_ other contractions. 


Here Ai are Lorentz scalars. In general, there are = (2n— 1 )!! different contractions. 
In many cases, the integral is symmetric under the interchange of some of its indices. 
Then one can symmetrize both sides of Eq. (2) such that the number of independent 
constants Aj is further reduced. Now, contracting 2 with each of the independent tensor 
structures in turn, one obtains a system of equations for the Aj with scalar integrals on 
the left hand site. Solving for the constants and reinserting the solution on the right 
hand side of 2, one obtains a decomposition of the integral in terms of scalar ones. 

The computation of the scalar integrals is described for each type of AI integral sepa¬ 
rately. Note however, that the Feynman rules for the composite operator 0{x) contain 
an additional factor of in comparison to the usual Feynman rules for an identical 
term in the Lagrangian (this results from the fact that there is no integration over x). 

The translation from the results of the Feynman diagrams to the operators occurring 
in the OPE can be done with different sophistication (for example with the help of 
the counterterm technique [17] or the method of projectors [18]). Here we simply give 
operators that yield vertices with the correct external fields and with vertex factors that 
are those polynomials resulting from the computation of the AI-graphs. 

The vertex factor corresponding to an operator 0(x) is: 

out in 

It is useful to note that the operator corresponds to the projection operator on the 
transversal component of the momentum that enters the vertex: 


5 _ 5 

SA-{p) “ SA-{p) 


= J d'^ke + C>(fls)| 

= - 5°-'’ {pagfj./3 - P/3flA.a) -t o (gs) ■ 


The generalization to higher powers of gluon fields is straightforward. 

Additional momentum factors are of course generated by additional derivatives of the 
gluon field operators. The situation is especially easy for additional factors of the overall 
momentum flowing into the vertex: 

n (4 m) («!:,.,<^)) 

= (-4+^6-*^"" |n +G(3.))|. 
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The C>(gs )-contribution is of higher order in comparison to the diagrams considered here 
and can be neglected. We shall denote the OPE as 

M—¥00 


O 


VC 

Hoop 


+ o. 


VC 

2loop 


+ O: 


VC 

3loop 


M- 


Q'y^Q 

Here indicates that the relation gives an asymptotic series. 


nTC ^TC ^TC 

^lloop ' ^2loop ' ^3loop' 


4.1. One-loop contribution 

4.1.1. Structure of the fermion loop 



Figure 2: Fermion loop with n gluon legs 

The graph 2 is the LO-contribution to the matrix element: 

Ga (*, XI, . . . , Xn) = (0| T O (x) (xi) . . . (Xn) |0) . 

Here 7“ (x) = (Q •y°‘TaQ) (x) is the color current. Now, in QCD the color current is not 
exactly conserved, but we have daJa (x) = 0 + G (gs)- The fermion loop that we consider 
contains only the contributions of LO in gs. Thus, to analyze the loop we can assume 
that the current is conserved. Taking the divergence of the matrix element, LO terms 
are only generated by moving the derivative inside the time ordering: 

dai Goj,'.;.“o„ (x, Xi, . . . , Xn) 

n n 

= ^ 5 (x? - x°) (0| |T O (x) [j°^ (xi) , Jal (xd] n j'j |0) 

i=2 

n 

+ S{x\- x°) (0| T [jo^ (xi) , O (x)] n ^7 ■ 

Now, the commutators are determined by the transformation properties of the operators 
under S'[/(3)c-transformations: 

[j°^ (xi) . G (x)] = 0 

[j°^ (Xl) , Ja- (Xi)] = ifa^aib (Xi) (xl - Xi) . 
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Here we assume that possible Schwinger-terms can be neglected for vacuum expectation 
values that we are computing. This gives the result 

n n 

da^Gal;::yn{xi,...,Xn) = - Xi) fa^aib {0\TO (x) (Xi) Ja^{xj)\0). 

i=2 j=2,ji^i 

In momentum space we have 

n n 

(Pl)ai {X,P 1 ,. . . ,Pn) = ^ fa^aib (0| T O (x) j“* (pi + p^) Ja^ (pj) |0) . 

i=2 j=2,j^i 

Obviously, there are similar relations for all external momenta P 2 • • -Pn- Now,we can draw 
conclusions for the case of three and four external gluons. 

For 3 external gluons, we have 

(Pl)ai (a:,Pl,P2,P3) = faia 2 b (0| TO (*) (pi + P 2 ) (Ps) |0) 

+ fa.asb (0| T O (X) + p 3 ) j“2 (^ 3 ) |o) 

= 0 . 


The last equality follows from C-invariance as was discussed above. From this relation 
(and from the corresponding ones for P 2 ,P 3 ) it is clear that 




.“a„ (a:,Pi,P2,P3) = n 




(x,Pl,P2,P3) ■ 


Now, the power counting arguments given above apply. These show that the matrix 
element for the tensor current is suppressed as M~^ and the vector current as M“^. 

For the case of four external gluons we have 


iPl)ai G'a(,...,an {x,Pl,P2,P3,Pi) 


fa.a^b (0| T O {X) J“2 (pi + P2) Ja| (ps) Ja^ (P4) |0) 


+ /aiagfc (0| T O {X) J“3 (pi + p 3 ) j“2 (^ 3 ) j“4 |o^ 

+ fa^a^b (0| T O {X) j;4 (pi + p 4 ) j“2 (^ 3 ) j“3 (^ 3 ) |o) 



Here n = 3 for the tensor current and n = 4 for the vector current. Since this applies 
for all external momenta, those parts of the matrix elements that are not suppressed as 
described above must have at least four projection operators. From power counting and 
the general structure of the expansion for tensor current and vector current we see that 
the contributions with these projectors are suppressed by M“^and M“®, respectively. 


4.1.2. Feynman integrals for the for the one-loop contribution 

The only 1-loop integrals to compute are shown in la and lb. The integrals were 
computed for a fixed assignment of external momenta. Afterwards the permutations in 
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external momenta were done to generate the crossed diagrams. For the diagram la we 
have the analytic expression: 


loop = - (igsf tro {Tann) I -0^ tv, (t 7° 






Here we use r = 7 ^ for the vector current and r = a^‘' for the tensor current. For the 
diagram lb with four external lines we have 


li\oop = - i {igst trc (TanTcTo) j -0 tr^ 7° 


$ l/i + 1/2 ~ ^ 


ft _ 1 _ r 1 p 1 1 

+ ^ + //[ + /^ — IS/I ^ ^ — M ^ ^ — M J 

The master integrals that occur after the expansion are: 


F(n) = 


= / 


(fk ( 1 




(_l)n. r(n-|) 

( 4 ^)d /2 r(n) • 


We use = —q — li — 12 to simplify the notation (this is the momentum entering in the 
last gluon vertex). 


4.1.3. Result for the vector current 

First consider the vector current in la. The result of expanding these diagrams up to 
the order l/M'^ reads: 

ilvc)i loop + perm. 

^ ImmC" ! I [ J 

+ {g^^q - g^'^q) {g^q - 5““ ih),) {gUg^ - g^^u.) {{h)^gl - sj (h),) 

+ {g^q - (^2).,) (<?£/ - ((h),flj - gl {Ii)s) 

-I [ - g'^^ij) (g^ig, - gW) (/'^2 - (sS {h)s - gf (h) J 

+ - g^’^q) (5“ (h).y - 5“ (h)p) {g^q2 - g^'^ii) {gUgs - g‘^g^) 

+ - g^’^q) [g^p ihq - g^ (h),) {g^^^ - g°^4) (g^gs - <?^-) ] }• 

This way to write the result clearly shows the expected projector structure of the result: 
Contracting with g,, {li)a, {h)p or {1^)t we obtain zero. 
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Due to the projector-structure of the results it is easy to write down the corresponding 
operators. After a short calculation, we have, omitting the argument x and neglecting 
higher terms in the coupling constant: 


ilvc)i loop + perm. 

S 5 S 

ih) ih) (Zi) 
S 5 S 

SA^ ih) SAb^ il2)SA<f (Zi) 


9s ^al a2 a3 


i-lGf,(x) (Go2)s, Gli^ + iGo2)s, G:1 


-(57 


72M^tv^ 


tTr 


Vp , !-{G^\ G^s}G^p + G^^Gs^GPP 
o 


With the application to instanton solutions in mind, we introduce the field strength 
tensor Fpo = gsGpy. Thus, the contribution can be generated by 


nvc 

^lloop 


-1 


72M47r2 


tVc 


V« 


,F^s]F^'' + Fs-^F^^PPP I + 


M6 y 


(3) 


This is exactly the result from the calculation of [5] after continuing to Minkowski¬ 
spacetime and rescaling the gluon fields. 

The contribution of figure lb will not be given here due to its huge size. However, it 
corresponds exactly to the O [plVcontribution to Eq. (3). 


4.1.4. Tensor current 


-igldqbc 

48M37r2 


For the tensor current a similar computation gives the result: 

{ItcYi loop + perm. 

[g^^ ih)p - g?^) - /D7) [g^r^ - g^'' ih)^) 

(<7“^ ii3)p - g?i^) - g’'^ij) {g^A - (^ 2 )^) 

{gP^lP - iPgP^) (g^Z^" - (Z2)^) (3“^ (Zg)^ - gpj) - ^ u 

- (ZZ“D^ - g^PFs) {iPg-^^ - g^^q) ( 3 ^ {h^ - 3^ (^ 2 ),) 

- - g^^q) - g^^q) {gp ih), - zz“ {h)p) 

- {qg’'^ - g^^ii) [g^qf^ - Z'’ (^ 2 ),) [gpl - 3 “"^ {h)p 

As expected the result vanishes for contractions with (Zi)a, {l 2 )g or (Za), 

Again, we can rewrite this in terms of an operator insertion: 


{ItcYi loop + perm. 

^_^_ [ -Z , 

SA%{h)SA<^{h)5A%{l2) "l24M3^2[ 


Ur Ur \j^p — Ur Ur \j^p Ur Ur-ypUr 
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Thus, we can write this contribution in terms of the operator insertion 


O 


TC 

Hoop 



-1 

24M3^ 


- F'^i^FP^'F^p + F'^PF-ypF^^'' 



( 4 ) 


which coincides with the result of [4] if one takes into account the analytic continuation 
to Minkowski-spacetime and the rescaling of the gluon fields. As in the case of the vector 
current, the diagram figure lb corresponds to the O [plVcontribution to equation 4. 


4.2. Two-loop contribution 


Since there are several different kinds of integrals at two-loop level, their explicit forms 
are given in the Appendices A and B. The integrals that have to be computed can be 
done in closed form and are of the type 


F{ni,n2) = J 


(_l)"l+"2+"3+i 


d‘^k d’^h _ 1 _ 

{2Tvf‘^ [fc2_M2]”l + hf - [i2]”3 

/ r (ni -I- 712 -t ns - d) r (n^ + ^3 - T (^2 + T 

y r (m-I-n2-I-2n3 - d) r (m) r (n2) r 



(5) 


As was discussed in section 3 there are no gauge invariant purely gluonic operators that 
can contribute to the order 1/M^ for the vector current or 1/M^ for the tensor current. 
Therefore, the l/M^(l/M)-contributions of the 2-loop diagrams should add up to 0 at each 
order in perturbation theory. At O and O the results given in the Appendices 
A and B confirm this reasoning; 

oXLp = o + 

O2 loop = 0 + 0 • 


4.3. Three-loop contribution 


In this case we have to compute the three loop integrals with external fermions (Ik). 
These are given by 


I3 loop — *(*<?«) [t'l'c (TaTfjTc) Ta Tf) Tc] J 


d‘^k d‘^li d^h f 1 1 


try r 


(27r)^^ I ^1 ^2 ( 9 + ^ 1 +^ 2 )^ 

/? 1 T 1 


’■-M ' fi + ^ + l/i+l^-M’ fi + ^ + l/i-M ' fi + fl-M 


^ -h ^ -h — m ^ + ^ + U + — ^ 
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Here m is the mass of any of the light quark flavors. 
The integrals that have to be computed are of the form 


/ {ai, 02, as, a4, a^, ae) 


r d'^kdHid% 



1 




1 {271^^ [fc 2 _ 

{k + hf - 1 

a2 

{k + h + hf -1 


{h + 

“6 


These integrals are hard to compute directly, but there is a way of computing them by 
the integration by parts (IBP) technique. Here we give a slightly modihed version of the 
approach of [19], where exactly these kinds of integrals where computed. One should 
also note that there are several software packages that compute such kinds of integrals 
(‘MATAD’ [20], ‘FIRE’ [21]). However, due to the fact that these integrals occur as 
an intermediate step after expanding in the heavy quark mass and tensor reduction of 
the result and that one has to do at the order of 50.000 of these integrals it was more 
convenient to write an own program that is specialized on exactly this type of integral. 


4.3.1. The IBP-technique 

In order to reduce the number of integrals to 
the symmetry properties of these integrals: 

/ (01,02,03,04,05,06) 

/ (01,02,03,04,05,05) 

/ (01,02,03,04,05,05) 

I (01,02,03,04,05,05) 

/ (01,02,03,04,05,05) 


compute, it is important to note some of 


= 7(02,03,01,05,05,04) 

= 7(01,03,02,05,05,04) 

= 7(02,01,03,04,05,05) 

= 7(03,02,01,05,04,05) 

= 7(03,01,02,05,04,05) 


These relations will be used such that the series of indices 01 , 02 , 03 is always decreasing 
unless one of the indices 04, 05,05 is non-positive and oi, 02, 03 are positive. In these 
cases the indices 04, 05,05 will be chosen to be decreasing. 

We will now show that we can use the IBP-technique to reduce an integral with arbitrary 
values of the indices o^ to integrals which fall in one of the following classes: 

1 . oi < 0 , 02 < 0 or 03 < 0 

2 . 05 < 0 and 05 < 0 . 

Afterwards, we show how to compute the integrals for these special cases. 
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Using the relations 


k ■ h 

1 

“ 2 



k-h 

1 

( 3 ") - i ( 2 “) + i ( 4 ") - i ( 6 “) 



2 

V / 2 \ J 2\J 2 \ J 


h ■ h 

1 

"" 2 




with the standard-notation 

I {■ ■ ■ ,ai,-■ ■) = I {... ,ai±l,... 


one can use 


0 = 


d’^kcThcrh d 


1^2 


( 27 r 


, 3 d 


[fe2 _ 1]-1 [(fc + _ 1] [(fc + + ^2)" - 1] ' 


[iir [{h+hfY 


to obtain 


I (ai,..., og) 


{o3 (3+) [(2-) - (5-)] +06 (6+) [(4-) - (5-)]} 
[as -I- 205 + ct^ — D] 


I{ai,..., ae) . 


This reduces 02 - 1 - 04 - 1 - 05 . The relation can be used until one of the indices is 0 . In the 
cases 04 < 0 or 05 < 0 this relation cannot be used. However, because of the conventional 
ordering of the indices, this case would imply a negative value of og and it would require 
05 to be the second non-positive index. If one finally arrives at 04 = 0 or 05 = 0 , one can 
stop, if og is already non-positive. If this is not the case, one can interchange the indices 
to make ag the index which is 0 . Then one can use the relation: 


7 (ai, 02 , 03 , 04 , 05 , 0 ) 


{03 (3+) [(2-)-(5-)]} 
[03 -I- 205 — 77] 


7 (oi, 02 , 03 , 04 , 05 , 0 ). 


This can be used until 02 or 05 is 0 . We conclude that in any case we end up with one 
of the special cases given above. 

Now we show how to compute the integrals for the special values of the coefficients. 


fTi < 0, 0-2 ^ 0 or 0^3 < 0! 

The cases a.\ < 0,02 <0 can be reduced to the case 03 < 0 . In this case one has the 
integral 


7 ( 01 , 02 , - I 031 , 04 , 05 , Og) 


d’^kd'^hd'^h 

(fc -1- hY -I- 27 

• I2 + 2^1 • I2 ^2 

Iq^sI 

{ 27 Tf‘^ [fc2 _ 1]«1 

\k + hf-l 


{h + hY 

ag 
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If one of the indices 05 or oq is nonnegative, the integral is 0 . Expanding the numerator 
and doing the obvious cancellations, one is left with integrals of the type: 


( f d‘^kd\ 1 

(/ 

d% 

{ 2 k-hf 

y ( 271 ) 2 "'+ 

l]“" [(2] “4 

[J 

( 271 )^ [(2] “5 

\ii+i 2 ?Yy j 


[§] 

= c ( 05 , 06 , 

r =0 


d^kd'^h 

)f] 

r 1 

r 91 ^^ 

(k + h] 

)2 - fc2 _ il 

p—2r 

( 271 ) 2 "' 

1]“ 


{k + hf - 1 

“2 |■^ 2 j“ 4 +“ 5 +"" 6-"'/2 


with 


j(_l)<i/2 r (^05 + 06 - I - r + I - 06 - r) r (^ - «5 + r) p! 
(47r)^/2 r (05) r (oe) r (p + d - 05 - Oe) r! (p - 2 r)! 


After the some cancellations, one has integrals of the type 

r d’^kdHi _ 1 _ 

J (2^)"^ [k 2 _ i]<-i +hf-i]"^ [iir^' 

If ai or 02 are non-positive, the integral is 0 . If all indices are positive, this is the 
standard-integral from the appendix. If 03 is negative, one can actually use the result 
( 5 ) as trivial analytic continuation to 03 < 0. 


05 < 0 and 06 < 0: 

Renaming and exchanging integration variables, we have: 


I (01,02, 03, 04, —05, —06) 


/ 


d'^kd\d'^l2 


(27 


\3d 


(TkdHi (Th 


(271 


,3d 



r, 2 i ^5 

'{h+hfY'' 

[k^-ifi y+h 


l“" [(fc+(i+(2)^-1]“^ [iir^ 

[^1 + ^2 + 2 Zi • (2! 

[^2 

-hfc 2 +2(2 • fc]“® 


[fc2 _ 1]“1 [^2 _ 1]“2 _ ;lJ«3 ^yo^A ' 


If any of the massive lines have non-positive indices, the integral is 0 . If 04 is non¬ 


positive, the integral reduces to a product of three integrals like / 
be easily done. If all indices are positive: 


2 L-k 


(277)"' [fc2-l] 


which can 


I (oi, 02,03,04) 

/fd‘^kd‘^li 1 /r d% [2/1 •(2]*’[2A:-(2]‘^\\ 

[J { 2 nf{ 2 .))' 

Doing the Z2-integration and simplifying one has the easy integrals (take p < q). 


I (oi, 02, O3, 04) 


p 

^D{aA,P,q,t) 

t=o 


d^^kdX [ll\'[k^]'^^'[2k-l,f 

( 271 )"^ [fe2-l]“l [(2-1] “2 + 
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Here 


D{a 4 ,p,q,t) 


I (P+9) 


(47r 


d/2 


r (a4 - (£+ 2 ) - I) 

f>H 


p!(p + g)! 

t\ip-2t)\(^+ty 


Note that in this case all terms in the numerator can be canceled since all exponents are 
integers. This gives; 

.. « ^ fd‘^kd‘^li 1 

I{a,p,^) = / - 575 - o -• 

J (27r)2‘^ 

This is again Eq. (5). 


4.3.2. Results for the three-loop contributions to the vector current 

The techniques described above could be used to compute the contributions for the 
vector current. We use os = |f to obtain: 

(Wc)3 *-p + Perm- = (v)'(f ” 2^(3)) 

( 2 I, u\ /33 10C(3)\\ 

I — - 11C(3) \-\q -<i q j > ( 16-3^ ) T 

Here we use antisymmetrization including the symmetry factors, for example: = 

\ ( 9 ^ 7 ^ — q'^1^)- In this case it is possible to read off the corresponding operators imme¬ 
diately. 

- md. - 11C(3)) + - d'^d'^) (T7.4>) (g - }. 

Note that we can freely interchange derivatives with covariant ones since the difference is 
of higher order in gs. Using the QCD equations of motion, y 4' = —im 4/ and T ^ = im T 
and the fact that the vector current for the light quarks is conserved, d^ = 0, we 

have 


^4'dT7^7p4'j = — m du ^ (4'7^'I') -I- O (gs) ■ 

This gives the result: 


O 


VC 
3 loop 



HM)! 
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4.3.3. Results for the three-loop contributions to the tensor current 


Here the computation yields the result: 

('rc)3,„„, + penn. = (fS)" _ C(3)) 

+ im(T^^ -3C(3)^ Q + C(3)^ j- 

In terms of operators, we have; 

- m^<7'^"M>(^^-3C(3)) + (i+ C(3)) |. 


One can use the equations of motion and standard relations for the 7-matrices to show: 


—V['9r7^7p5' = —-I-9^^ ^ ^'7^7 p' 5-t O (ps) . 

The result for the tensor current can now be written as; 


5. Result 


In the final result, we have to remember that each light flavor will contribute to the 
three-loop diagrams. Collecting the contributions, we have 


dp tr. [l{F^\F^s}F^P + Fsy'<F^P^ + O 


37 13C(3) 


48 


18 


727r2M4 


( 6 ) 


-(v)’5wS{--»/»/-'"'®/(|-V(3))+ jaiq*,,-'®/)} 

-rnkp''-' + F^F^F*”-) (7) 


Here denotes the sum over all flavors of quarks that are light compared to M. 

These operator identities can be used to analize intrinsic heavy quark content of light 
hadrons. One of immediate applications is the study of the intrinsic charm influence on 
the nucleon e.m. and tensor form factors. 
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6. Application: Heavy quark contribution to magnetic moment 
and electromagnetic radii of the nucleon 


We can use the results of the previous section to derive a relation for the gluonic con¬ 
tribution to the nucleons charge radius that is due to heavy quarks. We denote the 
momentum transfer by q = p' — p and the sum of the momenta by P = p' + p. First of 
all, we need the form factor decomposition for the vector and the tensor current. The 
decomposition for the vector current is standard: 

(p', cr'I ^7'"^(0) Ip, cr) = «(p', cr') -I-F 2 (^<J^^|w(p, cr) (8) 

The tensor current can be decomposed as [22]: 


(p^, (j'\ vi' Ip, (t) 


/ O f. ( 2\ , / 2\ , P^q'' -PU^ 

= u[p ,a ) Ht yq j ^ -- Et (q J 


2 m^ 


Ht (q^^^u{p,a) . (9) 


We actually need qi, (p',cr'| ^ |p,cr): 

qi, {p , a'\ ^ \p, a) 


= u [p , o') )^iq,^a^’'HT (^q^"j + 

= u (p', a') ^Ht (q^'j 


2 mjv 




2mjf 


Ht (q 


"N 

+ 


Here we introduced the tensor magnetic form factor kt + 2Hx and 

used the Gordon-identity, 


Additionally we introduce form factors of the gluon operator = tvc 


1{F^^^,F^s}E^'' + 


Fs^F^^F^^'' 
posed as t 


Since this operator is an antisymmetric tensor of rank 2, it can be decom- 
he tensor current. Thus, 


(p',(t'| s^''{q) Ip, a) 

= — im% u (j)', a') I ia^''S + 




2 mjv 


f 2 \ 

V / 


2mj^ 


•S (g^) U (p, cr) (10) 


and 


Qi' S^‘'{ 0 ) \p,a) 
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Now, we consider the nucleon matrix element of the vector current made of heavy 
quarks.Using the derived operator identity (6) we obtain the heavy quark contribution 
to the electric form factor and the magnetic form factor 

Gm = fQ (q^) + F^ (q^y. 



We see that if we neglect the chiral corrections of order ~ mf/mjM, than the leading 
~ contributions can completely be expressed in terms of singlet e.m. form 

factors. The subleading in heavy quark mass corrections of order ~ 1/M^ are expressed 
in terms of form factors of the gluon operator (10). 


This gives the following result for the heavy quark contribution to the electric and 


and 


\ 


magnetic radii, which we define as 

/Qs\3 5 


dG 


Q 


= - 6 - 


E/M 

— 


dq 




Q 


+ 


2 

rM 


-) 

6M2 

-Y 1 

‘2mj^ 

1 


247r2 M4 

as\3 

5 

-) 

6M2 

mf 


2m jq 


1 



37 13C(3) \ 

16 6 / 


2 16 


K^(0) -7/^(0) 




18 


Mat 


16 






2 ds 

dq^ 




( 11 ) 


(12) 
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Finally, the heavy quark contribution to the magnetic moment (in nucleon magnetons) 
is; 

(13) 

We see that the leading ~ 1 /M^ contribution to /ig has additional suppression by the light 
quark mass ^ mf/rriN indicating that the subleading correction ~ 1/M^ can be numerically 
more important than the leading one. The result for the heavy quark contribution to 
the magnetic moments (13) coincides with the result obtained in Ref. [6], the results for 
the e.m. radii (11,12) are new. 


From Eqs. (11,12), neglecting the corrections of order ~ mf/rriM, one obtains the following 
model independent limits for the e.m. radii of the nucleon: 


lim 


/ 2 ' 


/as ^ 


(37 

M—^oo 


'^Q 

\ TV z 

' 12 ' 

V 8 

lim 


/ 2 ' 

) = 

/as^ 


(37 

M—¥oo 


Q 

\ TT z 

' 36 ' 

^ 8 


13C(3) 

3 

13C(3) 

3 



A 

niNj ' 


(14) 

(15) 


We see that the limiting values of the heavy quark contribution to nucleon e.m. radii 
are always negative, so the intrinsic heavy quarks shrink the e.m. size of the nucleon. 
These limiting expression can be useful for lattice QCD simulations. 


Now we give numerical results for the heavy quark contributions to magnetic moment 
and e.m. radii in the case of the charm quark. We use the quark masses mu = 2.3 MeV, 
md = 4.8 MeV, = 95 MeV and me = 1.275 GeV [23] and the coupling constant as (m^'^ = 
0.3. The tensor form factors and the tensor anomalous magnetic moments were computed 
in the chiral quark soliton model [24, 25]. The values for the proton and neutron are 
HlpiO) = H^,J0) = 0.92, = H^ JO) = -0.27, = 4„(0) = 3.03 and 4^(0) = 

K^„(0) = 1.56 [26, 27]. Here isospin symmetry was assumed and the values were evolved 
from a renormalization scale of = 0.36GeV^ to = m^ by the method described in 
these publications. The numerical result for the charm quark contribution to the e.m. 
radii of the nucleon is: 



(^-5.1 + 54.8 (^i?(0)-S'(0) + 2S(0))) 10"®fm^ 



^-1.4 + 54.8 


^77(0) + 4m 


2 dS 



10"®fm^, 


(16) 

(17) 


Here we have the same result for protons and neutrons since the isospin violating effects 
due to mass difference of u and d quarks are negligible. We see that if the form factors 
of the gluon operator (10) are numerically of order unity, than the subleading term 
~ 1/M^ is larger than the leading term for the case of the charm quark. This stresses 
the importance of the estimates of the gluon form factors (10). One of possibilities to 
estimate these form factors is the theory of instanton vacuum [28, 29]. However, one 
can easily see that the gluon operator (10) is exactly zero on the instanton field. This 
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implies that one has to consider the contribution of the instanton-anti-instanon pair to 
this operator, hence the gluon form factors (10) are suppressed by the instanton packing 
fraction in the vacuum. Given such suppression, we can very roughly estimate the size 
of the gluon form factors as: 

gluon form factor ~ ~ 10~^. (18) 

-/Vc rt^ 

Here the factor (Nc - 2) takes into account that the gluon operator (10) is zero for the 
case of Nc = 2, the factor ~ ^ is the instanton packing fraction, and ~ 0.1 is the 
twist-4 contribution to the nucleon structure functions (see details in Refs. [4, 30, 31]). 

Surely, the estimate (18) is very rough. Nevertheless, it indicates that suppression of the 
gluon forma factors by the instanton packing fraction can be not enough to make the 
leading ~ 1/M^ term in charm quark mass expansion dominant numerically. Therefore, 
it is important to make more accurate calculations of the gluon form factors, we shall 
give detailed estimate of them elsewhere. Let us note obvious general properties of the 
gluon form factors (10): 1) they exactly zero for the case Nc = 2, 2) they are zero for 
general self-dual gluon field configurations. 

For the heavy quark contribution to the magnetic moment the leading ~ 1/M^ term is 
suppressed by the light quark mass ~ mf/mj^, therefore we expect that for the case of 
the charm quark the subleading contribution ~ 1/M^ is dominant numerically. Indeed, if 
one take the values of the parameters discussed above, one obtains the following result: 

HQ = (^ 2.0 • - 0.8 5 ( 0 )) 10 "^. 

Clearly the contribution of the gluon form factor is dominant, therefore it is very impor¬ 
tant to make an estimate of this form factor. 


7. Application: Heavy quark mass contribution to the tensor 
form factors 


The nucleon tensor form factors are defined by Eq. (9). The heavy quark contribution 
to the tensor form factors can be easily read of the operator identity (7) we derived. 


The corresponding equation contains the gluon operator ■— trc 




F-ypF^PF^'' 


, which nucleon matrix element we parametrize as: 


{pn'\ Ip, a) 

= m%u (p', a') I 


2im]\[ 




-t 




U [q ) >u{p,a) . 


(19) 
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The heavy quark mass contribution to the tensor form factors is less suppressed than 
corresponding contribution to the vector form factors - the leading contribution is of 
order ~ The result is: 




T > \ T, ) 36M 


as \3 5 


- - 4C(3)j nifE^ (<?2) + (i^/ {q^) + (<?^)) 




24M37r2 


V{q^) 


H^q^) = (—Y — 

^ ^ ^ y n j 36M 


4(<j2) = f^Y — 

’ y n j 36M 


- Y1 {q'^) - ( 9 ^) 

- (5 "^/^T («^) + (9^) 


24M37r2 


U{q^) 


‘‘N 


24M37r2 


U (q2) + 2t/ (g2) 


Neglecting the contributions suppressed by the light quark masses ~ ruj/mf^ we can 
derive the following heavy quark mass limits for the tensor form factors; 


)”.G) = (v)’ Gr ^ G) + («’)) + ° (Gf) 


M 


All these limits are expressed completely in terms of e.m. nucleon form factors, so these 
are model independent results. 


8. Conclusion and outlook 


Using method of expansion by subgraph we derived the heavy quark mass expansion of 
the vector and tensor currents up three loops. The corresponding operator identities are 
given by Eqs. (6,7). We applied these result to analyse the influence of the intrinsic 
charm on the nucleon e.m. and tensor form factors. We show that the leading orders 
of 1/M expansion are model independent and can be expressed in terms of known light 
flavour nucleon form factors. For the subleading terms one needs the calculation of 
the nucleon matrix elements of the gluon operators of dimension six. We showed that 
the corresponding operators in the instanton vacuum are suppressed by the instanton 
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packing fraction. However, still the contribution of the gluon operators can be important 
numerically for the case of charm quark. 

The heavy quark mass expansion Eqs. (6,7) can used to study the intrinsic heavy 
quark content of other hadrons, for example of vector mesons. The detailed studies of 
these topics will be published elsewhere. Also the methods developed here can be easily 
generalized to the case of other heavy quark made fermionic operators. For example, one 
can derive the heavy quark mass expansion of the twist-2 heavy quark quark operators. 
This gives us direct access to the intrinsic heavy quark content in the nucleon structure 
functions. 
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Appendix A. Two loop contributions at O ial) 


This appendix gives the results for the two loop integrals at O (section 4.2). 



Figure 3: Two loop contributions 


All diagrams were computed with a fixed assignment of external momenta. After the 
results were obtained, symmetrization in external momenta was done. Note that before 
symmetrization, the diagrams of Fig. 3b come with an additional factor of 1/2 due to 
the four-gluon-vertex. To simplify the expressions, we denote the three- and four-gluon- 
vertices as 




KL7 


= 9 ih - hr + ih - hr + 9 ^^ ih - wr } 

= - iff" ( 
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Then we have 


^2 loop 




(fk d'^r 


(27r 


iq + r + h,-q-li - l2,h - r) {-r,r - l2,h) 


2d aae 


beb 


tvc {TqTiTc) ~~n -n -9 

^ b J (g + r + /if [r-hf 


try ( r 


1 


^ + ^ + l/i_+f-M' V^ + ^ + li- M ' $ + M 

df'k df'r X 1 1 

b 


7 


1 


7 


^2 loop 


\ii9sf I 




{TaT-r-Tc) 4 - 2 

nhhn \ b / ^2 


1 


7 


1 


M ^ ^ + ^ + li+f-M' $ + ^ + l/i_- M ' ^ + 


I^iLp = iiWst j 


d<^k d^r ^rj,T (TaThTeT^) 


(27 


tr-y I r 


1 


\2d cec 


I(: + ^ +li+1/2-M' l(: + ^ + li-M 

1 r 1 


7 

1 1 
^ {r + hf 


7 


+ ^ + + / — Tf ^ ^ — A7 

- * / —-^ {TeTeTaTtTc) 

J 2nr‘^ 


jiv 

^2 /oop 


ir-y I r 


1 


i(:-M^^^ + f-M' $- M ' ^ + (^ + 1 / 1 + 1 / 2 -M 

1 r 1 

7 


7 


1 


7 


1 


We have three different integrals of type I, three different integrals of type II, six integrals 
of type III and ten different integrals of type IV. However, all integrals in one class, can 
in spite of there different topologies, be computed with the same methods. 


We computed all diagrams. For both the vector and the tensor current the sum of all 
contributions is 0. This is what was expected for reasons of gauge invariance. 
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Appendix B. Two loop contributions at 0{g!) 


This appendix gives the results for the two loop integrals at O (^ 5 ®^ (section 4.2). The 
computation here was very similar to that for the order gf. Thus, here only the different 
types of diagrams are given, not each diagram that was computed. 



(a) 2-loop diagram V 



(b) 2-loop diagram VI 



(d) 2-loop diagram VIII 



(c) 2-loop diagram VII 


Figure 4; Two loop contributions 


Again, the integrals were first done with a fixed assignment of external momenta. There 
are 6 diagrams of type shown in Fig. 4a, 10 diagrams of type 4b, 10 diagrams of type 
4c and 15 diagrams of type 4d. Finally, there are 24 permutations of external momenta. 
The explicit expressions for the diagrams are 


28 



Appendix C. Integrals necessary for the calculation of the 
three-loop integrals 

Here we give the results for the master integrals needed in Section 4.3.1. These integrals 
can be found in [32] and [33]. The notation |[gr]* with ‘2i + Yl'j=i ^ 

means the sum of all possible distributions of indices fii,... ,^n on i metric tensors and 
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momenta li. 


I (a,^ti. ..^n) 


I {a,P,p) 




r d‘^k fc^i ... fc^" 

J {2TTf [fc^ - 1]“ 

(47r)‘'/^ 2"/2 r (a) I ® i 

f d‘^k {2k ■ g)" 

J {27Tf ikT[ik+pff 

j(_l)d/2+n r(n+l) 

(4^)d/2 r(a)r(/ 3 )r(d + n-a-/?) 


H r (^a + /3r +n - a - r +i[2p • g]” 

^ r(t+ l)r(n- 2i + 1) [p2]«+/5-i-d/2 

r d^k d ‘^1 _ 1 _ 

J (27r)2‘^ _ 1]« |‘(fc + if _ ij ^ pf 

_(_^)a+/3+7 /r(a + /3 + 7 -d)r(a + 7-f)r(/3 + 7 -^)r(|- 7 ) 
(47r)‘^ r(a + ;3 + 27 -d)r(a)r(/ 3 )r(^) 
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